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1. Introduction 

In recent years many studies have been devoted to nonequilibrium statistical-mechanical 
models yielding condensation, such as urn models [1, 2, 3, 4], zero-range processes 
(ZRP) [5, 6, 7, 8, 9], generalised mass transport models [10], or driven diffusive 
systems [11, 12, 13, 14]. In all these models except for the latter, the condensate 
manifests itself by the macroscopic occupation of a single site of a thermodynamically 
large system by a finite fraction of the whole mass. For driven diffusive systems, the 
condensate manifests itself as a domain of macroscopic size. 

For a large but finite system at stationarity, the condensate is expected to perform 
some ergodic motion over the sample. The aim of this work is to investigate in detail 
the nature of this motion in the stationary state, on the example of a class of ZRP. We 
demonstrate that this motion is fully non-local, and we determine the corresponding 
characteristic time r. In the situation of most current interest where the ZRP has a 
power-law critical mass distribution, we find that r scales as a power of the system 
size, in any dimension. Finally, our predictions are generalised to other types of critical 
mass distributions. Whatever this distribution, r never scales exponentially with the 
system size. Little attention had been devoted so far to the stationary-state motion of 
the condensate [8, 10]. A critical reading of these references is given at the end of the 
present work. 

2. Definition of the model 

A ZRP is a stochastic process, defined as follows [15, 16, 9]. Consider a finite connected 
graph, made of M sites. Particles hop from site to site with rates which only depend 
on the occupation of the departure site. Denote by iVj the instantaneous occupation of 
site i (i — 1, . . . , M). The total number of particles Ni + ■ ■ ■ + N M = N is conserved by 
the dynamics. 

An elementary step of the dynamics consists in choosing a departure site d and an 
arrival site a connected to site d, and in transferring one of the particles present on site d 
to site a. This process takes place with rate Uk per unit time, where k = Nj ^ 0. In the 
mean-field geometry, all sites are connected, i.e., sites d and a are chosen independently 
at random. On finite-dimensional lattices, site a is chosen among the first neighbours 
of site d, possibly in an asymmetric way. In one dimension, the most general dynamics 
corresponds to choosing site a to be the right neighbour of site d with probability p, 
or its left neighbour with the complementary probability q — 1 — p. In the following 
we only consider the one-dimensional symmetric dynamics, corresponding to p — 1/2, 
and the fully asymmetric one, corresponding to p — 1, both with periodic boundary 
conditions. 

The fundamental property of ZRP is that the stationary distribution of particles 
is factorised, and explicitly known in terms of the rates [15, 9]. The probability 
of a given configuration of the system {iVi, . . . , N M } is equal to a product of weights 
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associated with the occupation of each site (see Appendix A for more details): 

l M 

«) = — 11m, (2-1) 

where 



Z M ,N i=1 



Po = 1, Pk = • (2.2) 

Ui...U k 

The normalisation factor Zm,n plays the role of a partition function. The stationary 
occupation probabilities of any given site read 

„ PkZM-l,N-k / qA 

Jfc.st - ^ • (^.O) 

In the present work, unless otherwise stated, we consider the rate [6, 11, 7, 8] 

u k = 1 + p (2.4) 

where b is a given parameter, which plays the role of inverse temperature. With this 
choice of rate, in the thermodynamic limit (M — > oo at fixed density N/M = p), the 
system has a continuous phase transition at the critical density 

1 

whenever b > 2, i.e., at low enough temperature. The critical density separates a fluid 
phase from a condensed phase. In the fluid phase (p < p c ) , the occupation probabilities 
fk st fall off exponentially. At the critical density (p = p c ), they fall off as a power law: 

A,, « P-5) 

(see (A.8)). In the condensed phase (p > p c ), for a large and finite system, the 
particles form a uniform critical background and a macroscopic condensate, consisting 
(on average) of A excess particles with respect to the critical state, where 

A = N-M Pc = M(p- Pc ). (2.6) 

The condensate appears as a hump in the stationary distribution /fc ;St (see Figure 3 
below). 

The case of the rate [11, 12] 
a 

where a is an arbitrary exponent, is discussed at the end of the present work. 



3. Heuristic approach 

For a large and finite system in the stationary state, as time passes, the condensate keeps 
moving across the system: it spends long lengths of time on a given site, before suddenly 
disappearing and reappearing on another site. The typical value of these lengths of time 
defines the characteristic time r of the dynamics of the condensate. The aim of this 
work is to characterise how r scales with the system size M. 
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An intuitive understanding of the phenomenon is easily gained by performing simple 
Monte-Carlo simulations. Such simulations done in the three geometries introduced 
above: mean-field (MF), one-dimensional asymmetric (IDAS) (p = 1), and one- 
dimensional symmetric (IDS) (p — q — 1/2), lead to a common picture. 




3 4 

t [xlO 4 



Figure 1. Dynamics of the condensate (IDS geometry with b = 4, M = 40, 
N = 80). Upper panel: instantaneous number of particles iV max (i) on the most 
populated site. Lower panel: location i m &x(t) of that site. 



The condensate is immobile for rather long lapses of time; it then performs sudden 
random non-local jumps all over the system, at Poissonian times whose characteristic 
scale grows rapidly with the system size M. Figure 1 illustrates this process for the IDS 
case, for a system of size M = 40, with N = 80 particles, i.e., p — 2, and 6 = 4, hence 
p c = 1/2. The upper panel shows the track of the instantaneous number of particles 
N max (t) on the most populated site. The signal for N max (t) fluctuates around A = 60, 
the mean size of the condensate (see (2.6)). The lower panel shows the label i max (t) of 
that site, i.e., the location of the condensate. The non-local character of the motion of 
the condensate is clearly visible, whereas the longest lapses of time where the condensate 
stays still give a heuristic measure of the characteristic time r. 

These features will be turned to quantitative predictions in the following. We 
first analyse in detail the case of two sites (M = 2) (Section 4). We then introduce an 
approximate model for the dynamics of the ZRP, for M arbitrary, based on a Markovian 
Ansatz (Section 5). The relaxation time r MA within the Markovian Ansatz will be shown 
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to grow asymptotically according to (5.5). Its dependence on M and p is the following: 

A b+1 

T^~{p- Pc ) b+1 M b ~—. 

We finally present the results of Monte-Carlo simulations for the temporal correlations 

of local and global density fluctuations in the stationary state of the ZRP (Section 6). 

The associated respective relaxation times T\ oc and r g i will be shown to scale 

the MF and IDAS geometries, and as Mr MA ~ A 6+1 in the IDS geometry, up to finite 

prefactors. 

For clarity we have regrouped in two appendices the relevant tools for the derivation 
of the results of sections 4 and 5. Appendix A is concerned with the static properties of 
ZRP, Appendix B with the analysis of the characteristic times appearing in the biased 
motion of a random walker on an interval. 



4. The case of two sites 



The case of two sites (M = 2, N arbitrary) is exemplary because it shares a number 
of features with the general case. In the stationary state, for N large, the system flips 
between two states where most of the particles are either on one site, or on the other one. 
In this section we determine the flipping time of the system, that is the characteristic 
time for the system to flip between these two states. 

The dynamics of the two-site system is encoded in the master equation for the 
occupation probability of site 1, defined as fk(t) = V(Ni(t) = k). The temporal 
evolution of this quantity is given by the equations 
df k (t) 



dt 
d/o(t) 

dt 
df N (t) 



= Hk+i fk+i + f k -i - (p k + \ k )f k (1 < k < N - 1), 
= Afi/i-A / , (4.1) 
— Ajv-i /at-i — ^nJn, 



dt 

where = UN-k and pk = Uk are respectively the rate at which a particle enters site 1, 
coming from site 2, or leaves site 1 for site 2. The equations for k = or k = N 
are special, since uq = 0. The above equations describe a biased random walk on the 
interval (0, N), with reflecting boundaries at and N, the position of the walker being 
the random variable N\(t), i.e., the number of particles on site 1. 

We begin by a description of the statics. The time-independent solution to (4.1) is 
easily seen to coincide with the prediction of the general formula (2.3) for the stationary 
probabilities 

, _ PkPN-k ( , 9 , 

Jk,st — , 

where pk is given by (2.2) and Z 2) n by (A. 3). Using the asymptotic form (A. 6) of p k at 
large k, the probabilities (4.2) are found to scale as 

/M ~ k b( N _ k y 
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for large k and N — k. The amplitude Cn is given below (see (4.3) and (4.4)). The 
distribution profile is therefore U-shaped for any positive b: the most probable stationary 
configurations are those where almost all the particles are located on one of the two 
sites. In the language of the effective potential, defined as — ln(/fc )St //o,st), these two 
states correspond to two symmetric valleys, on both sides of a potential barrier situated 
at N/2. However, as we now discuss, the cases b > 1 and < b < 1 are qualitatively 
different. They correspond respectively to the case of a high and a low potential barrier 
(see Appendix B). 

For b > 1, the partition function Z 2j n is dominated by configurations such that 
almost all the particles lie on either of the sites, i.e., either k is finite, or N — k is finite. 
Indeed the sum of the p k , i.e., -P(l) (see (A. 7)), is convergent. We have 

^^^m+v^, c^ftim*. ,4.3) 

For < b < 1, on the contrary, the sum of the pt is divergent. As a consequence, 
all the values of k contribute to the expression (A. 3) for the partition function Z2,n- 
Evaluating this expression by means of an integral, we obtain 

^ r(fe + iffli-fe) 2 r(2-2&) U4) 

A natural definition of the flipping time is the crossing time T = T L = T R , i.e., the 
mean time taken by the walker to cross the system from right to left, or from left to 
right. As shown in Appendix B, the expression of T only relies on the knowledge of the 
stationary state. It is given in terms of /fc jBt by (B.8), from which analytical asymptotic 
estimates of T are easily obtained. Replacing the sum by an integral in (B.8), we have, 
for b > 1 (regime of high barrier) , 

&r(6 + i) , +1 

T ~ (b-l)T(2b + 2) N ' (45) 
while for < b < 1 (regime of low barrier) , 



7rb 



N 2 . 



2(1 -46 2 )tan7r6 

We complement this study by the determination of the spectral properties of the 
master equation (4.1). The latter reads formally 

^ = EMkMt), (4.6) 

where Ai is the Markov matrix of the process. Denoting its eigenvalues by E a and the 
corresponding eigenvectors by <fk,a, the eigenvalue problem for Ai consists in solving 
the equation 

Ufc+lVfe+l,o + WjV+l-feVfc-l,o - (Uk + UN-k)<-Pk,a = ~E a Lp k ,a, (4.7) 

obtained by looking for solutions to (4.1) of the form fk{t) = ^pk^~ Et - The ground state 
£"0 = corresponds to the stationary solution tp^Q = /fc st of (4.2). The inverse of the 
other eigenvalues, r a — 1/E a for a — 1, . . . , N, are the characteristic times of the model. 
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Let us analyse the scaling behaviour of these times for N large. We make use of 
the results of Appendix B for the case where the potential has two symmetric valleys. 
For b > 1, the main characteristic time T\ is related to the crossing time T by (B.14): 



T 



(4.8) 



It thus scales as N b+1 , with an explicitly known prefactor (see (4.5)). The subsequent 
characteristic times r 2 , r 3 , . . ., diverge only as the diffusive timescale N 2 , i.e., much less 
rapidly than t± for b > 1. This can be shown by taking the continuum limit of (4.7). 
Introducing the scaling variable x = k/N and £ = N 2 E, we get the following differential 
equation 

x 2 (l-x)V" +bx(l-x)(l-2x)ip' + [£x 2 (l-x) 2 -b(l-2x + 2x 2 )]<p = 0(4.9) 

(a confluent Heun equation) for the scaling function tp(x), with Dirichlet boundary 
conditions: ip(0) = ip(l) = 0. Let £ a be the eigenvalues (a = 2, 3, . . .). We have 

N 2 

Ta ~ (4.10) 

A numerical integration of (4.9), with 6 = 2, yields £ 2 = 26.815, £ 3 = 74.609, and so on. 




Figure 2. Log- log plot of the crossing time T and of the first three characteristic 
times ti, T2, T3, against the number iV of particles, for the two-site problem 
with 6 = 2. Full straight lines: asymptotic scaling laws (4.5), (4.8), and (4.10), 
including prefactors: T iV 3 /30, n w iV 3 /60, r 2 w A^ 2 /26.185, r 3 w 
iV 2 /74.609. 



Figure 2 shows a log-log plot of the exact crossing time T, as given by (B.7), and 
of the first three characteristic times t±, r 2 , r 3 , obtained from solving (4.7) numerically, 
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against the number N of particles, for 6 = 2. The agreement with the analytical 
asymptotic estimates (4.5), (4.8), and (4.10) is excellent. 

For < b < 1, all the characteristic times, Ti, T2, . . ., grow proportionally to the 
diffusive timescale N 2 , as in (4.10). 

Consider finally the two-time correlation function for the occupation, computed in 
the stationary state: 

C(t) = (N 1 (t)N 1 (0))-(N 1 ) 2 , 

with (Ni) = N/2. This correlation function is the simplest example of an observable 
falling off to zero as C(t) ~ exp(— t/ri), where T\ = 1/Ei is the main characteristic 
time. Indeed, in the present case, the temporal evolution of C(t) follows the Markovian 
scheme presented above. We have 

(N^t) jVx(O)) = Yl JkViN^t) = fc|iVi(0) = j) V{NM = j), 

j,k>l 

i.e., 

(N^t) Aq(0)> = £ j k Vhk (t) f hSt = k lk (t), 

j,k>l k>l 

where the quantities 

p jtk (t) = ViN^t) = fe|JV x (0) = j), 7fc (t) = J Pjdt) kst 

follow the master equation (4.1), with initial conditions 

Pj,k(0) = S jtk , 7fe(0) = fc/fc,st, 

consistently with C(0) = Efe ^ 2 /fe,st - (A^i) 2 . 

To summarise, the case of two sites is instructive in several respects. The model 
exhibits two qualitatively different dynamical regimes. For b > 1, the dynamics is an 
example of a diffusion in a potential with a high barrier. The flipping time, namely the 
crossing time T, diverges as N b+1 , as does the main characteristic time r±. The limit 
value of the ratio T jr x is equal to 2, characteristic of a symmetric potential with a high 
barrier. For < b < 1, the flipping time, as well as all characteristic times, now grow as 
the diffusive timescale N 2 . The ratio T/t\ is non trivial and depends on b. Appendix B 
treats the general case of a non-symmetric potential, for use in the next section. 



5. The general case: Markovian Ansatz 

We now proceed to the analysis of the stationary dynamics of the condensate for the 
general case (M and N arbitrary, b > 2, p > p c ). We first make a crucial observation on 
the behaviour of the occupation probabilities f k , s t m the region k ^> 1 and A — k ^> 1, 
hereafter referred to as the region of the probability "dip". We refer the reader to [10, 9] 
for a general study of the statics of mass transport models with critical background, for 
a finite system, as well as to the additional information contained in Appendix A. 
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Figure 3. Logarithmic plot of the occupation probabilities fk >st in the 
condensed phase (b = 4, p = 4p c = 2), against the ratio k/M. Top to bottom: 
M = 20, 40, 80, 160, and 320. Full lines: / fcjSt obtained by (A.l) and (A.2). Full 
(empty) symbols: maxima (minima) of occupation probability. Dashed vertical 
lines: asymptotic locations of the minima: k/M = A/(2M) = {p — p c )/2 = 3/4, 
and of the maxima: k/M = A/M = p — p c = 3/2. 



Figure 3 shows a logarithmic plot of the occupation probabilities f^ s t 111 the 
condensed phase, computed using equations (A.l) and (A.2), against the ratio k/M, 
for b = 4, p = 2, and several values of M. This plot exhibits the following features. For 
k <C M, the distribution fk )S t is approximately given by the power law (2.5) of an infinite 
critical system. The contribution of the condensate slowly builds up as a probability 
hump around k — A, the mean number of excess particles. Finally one observes a broad 
and shallow probability dip in the region located between the critical background and 
the condensate hump. It is shown in Appendix A that the region of the dip is dominated 
by configurations where the excess particles are shared by two sites (see (A. 11)): 

The observed locations of the maxima (k « A) and minima (k « A/2) of the occupation 
probabilities corroborate this picture, as explained in the caption of Figure 3. 

These observations lead us to a first crude estimate of the characteristic time of the 
stationary-state dynamics of the condensate. The corresponding decay rate can indeed 
be argued to be proportional to the minimum / min of the stationary probabilities fk jSt . 
We thus obtain 

1 , , 
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We shall refer to this kind of estimate as the Arrhenius law, as it generalises the formula 
originally due to Arrhenius [17] for a chemical reaction rate in terms of the corresponding 
activation energy. In the present context, the result (5.1) implies that / m i n is reached 
near the middle of the dip region (k ~ N/2), and that the Arrhenius law (5.2) yields 

r A ~ A 6 . (5.3) 

We now turn to the derivation of a more refined approximation, yielding the correct 
prefactor to the Arrhenius law, somewhat along the lines of the approach initiated by 
Kramers [18] for chemical reactions. The above observations on the statics suggest to 
propose an effective description of the stationary dynamics of the condensate, based on 
a Markovian Ansatz. Assume that the condensate is on site number 1 at the initial 
observation time (t = 0). The number Ni(t) of particles on that site is initially very 
large, iVi(O) rs A, and therefore evolves slowly, until the condensate dissolves into the 
critical background. Thus 

• We single out Ni(t) as the collective co-ordinate of the system, that is the 
appropriate slow variable describing the dynamics of the condensate. 

• We model the dynamics of Ni(t) by (B.l), i.e., by a biased diffusive motion on the 
interval k = 0, . . . , N. The left hopping rates are taken equal to the microscopic 
ones: = The right hopping rates are so chosen that, in the stationary 
state, the probabilities fk, s t of the effective model coincide with the occupation 
probabilities (2.3) of the original ZRP. The detailed balance condition (B.2) yields 

^ _ /^fc+l/fc+l.st _ ZM-l,N-k-l 
fk,st ZM-l,N-k 

The rate is thus a function of k, M, and N. For M — 2, this formula gives 
Afe = UN-k, as expected. In the fluid phase, in the thermodynamic limit, the 
rates converge to zo, defined in (A. 5). Finally, in the condensed phase, in the 
dip region, we obtain 

Afe " 1 + A^" 

This effective description therefore reduces the full model to a Markovian model 
for one degree of freedom in an asymmetric potential. The two valleys of the potential 
are separated by a high (power-law) barrier. The left potential valley, corresponding 
to the critical background, has a weight Pl ~ 1, whereas the right potential valley, 
corresponding to the hump of the condensate, has a weight Pr ~ 1/M <g; 1. 

In this framework (see Appendix B) the stationary dynamics of the condensate 
is characterised by a single diverging timescale. We choose to define this timescale, 
denoted by r MA , to be the crossing time T L from the right valley to the left one in the 
effective Markovian problem. The characteristic time is thus expressed by (B.7), i.e., 

N j N 

T MA = T L =Y^ 7 fm,st, (5-4) 

e=i Me, st m= i 

in terms of known quantities, the rates and the stationary probabilities /fc )St . 
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Its asymptotic growth is easily determined by noting that (5.4) is dominated by 
the behaviour of the probability fk tS t in the region of the dip. Hence, inserting the 
expression (5.1) into (5.4), and evaluating the sum as an integral, we obtain 

6r(6 + l) A b+1 6r(6 + l) 



7"ma 



(P - Pc) b+1 M 



(5.5) 



(6-l)r(26 + 2) M (6-l)r(2& + 2) 
This result incorporates the full prefactor to the Arrhenius estimate (5.3), including a 
factor of A/M — p — p c , as well as an absolute numerical factor. It is also worth noticing 
the exact correspondence between (5.5) and the formula (4.5) for the two-site system, 
up to the replacement of M by 2 and of A by N, and up to an absolute factor 2. 





Figure 4. Log-log plot of the prediction r MA of the Markovian Ansatz for the 
characteristic time of condensate dynamics, against the system size M, for 
6 = 4 (hence p c = 1/2) and several values of the density in the condensed 
phase. Symbols: exact values of r MA , obtained by means of (5.4). Full straight 
lines: asymptotic scaling result (5.5), including prefactor. Dashed straight line: 
diffusive timescale of the critical dynamics. 



Figure 4 shows a log-log plot of the characteristic time r MA , given by (5.4), against 
the system size M, for b = 4 and several values of the density in the condensed 
phase. The scaling form (5.5) only sets in for rather large systems. For p — p c 
small, there is indeed a crossover phenomenon between the diffusive scaling of the 
timescale at the critical point, r MA ~ M 2 , irrespective of b, and that of the condensed 
phase, t ma ~ (p — p c ) b+1 M b . The crossover takes place for a system size of order 
M* ~ (p-p c )-( 6+1 )/( 6 - 2 ). The apparent crossover size M+, beyond which the form (5.5) 
holds, remains however rather large, even when the density p largely exceeds p c . For 
instance, for b = 4 and p = 4p c = 2 (the values used in the Monte-Carlo simulations of 
the next section), one has M* m 40. 
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It is also worth considering T R , the other crossing time of the effective Markovian 
problem. First of all, the estimate (B.10) for Tr fails in the present situation. Indeed, 
this estimate is derived in the framework considered in Appendix B, where the effective 
potential is assumed to be small everywhere except in the barrier. The present form 
of the right valley is however very different. The stationary probabilities fk, s t indeed 
fall off very fast as k goes to its maximal value N, so that the effective potential V k 
exhibits a sudden rise as k — > N . As a consequence, for a large system size, the exact 
formula (B.7) for T R is not dominated by the region of the dip. It is instead entirely 
dominated by the last stationary probability, which reads /jv, s t = Pn/Zm,n- Using the 
expression (A. 10) for Zm,n, we get 



The crossing time T R therefore grows exponentially with the system size M. The 
corresponding activation energy reads — \nf NtSt — V N « M\n(b/(b — 1)) = —M!F C , 
where T c is the critical free energy density. This crossing time has a clear physical 
interpretation: it is the characteristic time of the occurrence of the very improbable 
event where the whole critical background of the ZRP collapses onto the condensate 
(whence the occurrence of J- c ). The occupation of the condensate is then equal to 
N, well above the usual A. The exponentially large timescale T R is absent from the 
spectrum of the Markov matrix of the effective description, whose first eigenvalue scales 
as Ei ps 1/(P L T L ) ~ l/T L . The occurrence of exponentially large first-passage times 
has been underlined in several instances of urn models [4, 19]. 

6. The general case: numerical simulations 

In order to test the predictions of the Markovian Ansatz, we now compare them to 
results of numerical simulations. We have performed Monte-Carlo simulations of the 
ZRP with the rates (2.4), using random sequential updates, in the three geometries 
defined above: mean-field (MF), one-dimensional symmetric (IDS) (p — q — 1/2), and 
one-dimensional asymmetric (IDAS) (p = 1). We have focussed our attention onto the 
following stationary-state correlations, which are devised to be especially sensitive to 
the dynamics of the condensate: 

• Local density correlation function: 



where (Ni) = p. The initial value Ci oc (0) = Y,k=o k 2 fk,st — P 2 reads approximately 




CUt) = (^i(0^i(0)) - (iVx) 2 , 




Cgi(*) = i^*(0)>, 
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where 

N 

g(t)=J2N m (t)e 2 ^ M 

m=l 

is the first Fourier coefficient of the density fluctuations. This definition only makes 
sense for the one-dimensional geometry with periodic boundary conditions. 

Both correlation functions fall off to zero in the limit of long times, by construction. 
The correlation functions measured in simulations are found to obey a very clear 
exponential fall-off, except for very small systems, confirming thus the emergence of a 
single well-defined slow timescale characterising the dynamics of the condensate. In the 
one-dimensional geometries (IDS and IDAS), the local and global density correlations 
have different relaxation times: 

Q oc (t) ~ exp(-t/r loc ), C g i(t) ~ exp(-t/r g i). 

In the mean-field geometry only t\ oc can be measured. 

12 I 1 i 1 i 1 i 1 i 1 7~i 



o 
o 
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In M 

Figure 5. Log- log plot of the relaxation time T\ oc measured from the fall-off 
of the local density correlation Ci oc (i), against the system size M (b = 4, 
p = 4p c = 2). Symbols: numerical data for the three geometries considered. 
Full line: prediction t M a of the Markovian Ansatz. 



Figure 5 demonstrates that the prediction of the Markovian Ansatz correctly 
describes the main dependence of the measured relaxation time t\ oc of local density 
fluctuations (the relaxation time r g i has been omitted for clarity). The agreement is 
especially visible in the MF and IDAS models, whereas the relaxation time of the IDS 
model grows slightly faster. The prediction r MA of the Markovian Ansatz therefore 
provides a good estimate of the characteristic time of the dynamics of the condensate, 
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even for moderately large system sizes, below the crossover size M* where the power-law 
scaling result (5.5) sets in. 

In order to analyse in more detail the behaviour of the relaxation times t\ oc and r g i 
of local and global density fluctuations, we introduce the dimensionless time ratios 

-Rioc = — — , -Rgi = — — ■ (6.1) 

7"ma Tma 



2.0 




0.5 - 
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M 

Figure 6. Plot of the relaxation time ratios R\ oc and R g \ introduced in (6.1), 
against the system size M, for the MF and IDAS geometries (b = 4, p = 4p c = 
2). Symbols with error bars: numerical data. Full lines: Rational fits (see text) 
leading to the limiting values « 2 (MF) and E^} » R { ™ ] » 1.26 (IDAS). 



Figure 6 shows a plot of the time ratio R\ oc in the MF model, and of both time 
ratios i?i oc and R g \ in the IDAS model. The three series of data remain of order unity, 
but exhibit a weak residual growth as a function of M, at least in the range of system 
sizes shown on the Figure, which slows down beyond system sizes of the order of the 
crossover size M* 40, and seems to eventually saturate. Let us stress the difficulty in 
measuring these time ratios with a statistical error of a few percent for larger systems, 
whereas the measured relaxation times vary over four orders of magnitude. The fact 
that the time ratios saturate to finite values in the thermodynamic limit is further 
supported by the fits of the data by fractional linear functions of 1/M 2 , shown as full 
lines on Figure 6. The rationale behind this choice of a fitting function is the following. 
For 6 = 4, r MA scales as M 4 , whereas subleading terms proportional to the diffusive 
timescale M 2 are expected, among other correction terms. The presented fits lead to 
the following limiting values of the time ratios: ~ 2 (MF) and ~ R^ ~ 1.26 
(IDAS), hence 

7~loc ~ 2 T MA (MF), ri oc «r g i«1.26r MA (IDAS). (6.2) 
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The scaling formula (5.5) of the Markovian Ansatz prediction therefore correctly 
describes the asymptotic behaviour of the relaxation times in the MF and IDAS cases, 
up to finite prefactors. 



12.5 




Figure 7. Plot of the relaxation time ratios R\ oc and R g \ introduced in (6.1), 
against the system size M, for the IDS geometry (b = 4, p = 4p c = 2). Symbols 
with error bars: numerical data. Full straight lines: least-square fits with slopes 
0.12 and 0.30. 

Figure 7 shows plots of the time ratios R\ oc and R g \ in the IDS geometry. Both ratios 
exhibit a clear linear dependence on the system size M, implying that the relaxation 
times of local and global density fluctuations scale as 

n oc « 0.12 Mr MA , Tgi ~ 0.30 Mt ma (IDAS). (6.3) 

These expressions of the relaxation times in the IDS geometry therefore contain an extra 
factor of the system size M, with respect to the analytical prediction of the Markovian 
Ansatz, and to the results in the MF and IDAS geometries. This factor originates in 
the Gambler's ruin problem [20] (see Discussion). To close up this part on numerical 
results, let us emphasise that the prefactors entering (6.2) and (6.3) have been extracted 
accurately only for 6 = 4 and p = 2. These prefactors are however expected to have a 
weak dependence on b and p. 

7. Discussion 

This paper is devoted to the dynamics in the stationarity state of the macroscopic 
condensate occurring in a class of nonequilibrium statistical mechanical models. The 
problem is well-defined for a finite system consisting of M sites, where the condensate is 
made of a large but finite number A = M(p — p c ) of excess particles, up to fluctuations. 
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In this work, we have introduced, for the first time, a description of the stationary- 
state dynamics of the condensate, on the example of a class of ZRP. This effective 
description, based on a Markovian Ansatz, relies deeply on the exact knowledge of the 
stationary- state occupation probabilities /fc,st of the finite system. The essential feature 
of these probabilities is the occurrence of a dip region dominated by configurations 
where the excess particles are shared by two sites. The effective Markovian problem 
thus derived is an example of biased diffusion in the presence of a bistable potential 
with a high barrier separating two valleys. The present approach turns the Arrhenius 
estimate r A ~ l// m in of the characteristic time into the accurate prediction r MA , just 
as the method initiated by Kramers improves the Arrhenius law in the context of 
chemical reactions. Our main result is the following. The characteristic timescale for 
the dynamics of the condensate, r MA , always grows much faster than the diffusive time- 
scale M 2 , but not exponentially fast. For the class of ZRP with static exponent b > 2, 
the characteristic time r MA grows as M b . 

These theoretical predictions are corroborated by the results of Monte-Carlo 
simulations in the mean-field (MF), one-dimensional asymmetric (IDAS), and one- 
dimensional symmetric (IDS) geometries. Local and global temporal correlations of 
density fluctuations in the stationary state have been respectively measured by the 
correlation functions Ci oc (t) and C g \(t). The corresponding relaxation times ri oc and r g i 
are found to scale as follows: T\ oc ~ r MA ~ M b (MF), ri oc ~ r g i ~ r MA ~ M b (IDAS), and 
r ioc ~ r gi ~ Mr MA ~ M b+1 (IDS), up to absolute prefactors of order unity, which are 
expected to bear a weak dependence on b and p. The local and global relaxation times 
seem to be asymptotically equal (T g i/Vi oc « 1) in the IDAS geometry, and asymptotically 
proportional (t & \/t\ oc 2.5) in the IDS geometry. This observation confirms the visual 
impression gained from Figure 1: the stationary dynamics is non-local, in the sense that 
the jumps of the condensate spread more or less uniformly over the whole sample. (Any 
more local kind of dynamics, where large jumps are either absent or improbable, would 
imply that the ratio t^\/t\ oc diverges with the system size.) 

It is of interest to put our results on the dynamics of the condensate in the stationary 
state in perspective with those of the non- stationary dynamics, that we first recall. The 
non-stationary dynamics of the formation of the condensate has been investigated in 
some detail for ZRP, both in the mean-field geometry [7] and in one dimension [7, 8]. 
The following scenario holds in the MF and IDAS geometries. In the first stage of the 
dynamics, a coarsening phenomenon takes place: the typical occupation of the most 
populated sites grows as t 1 / 2 , whereas the number of those sites decays as M/t 1 / 2 . After 
a time of order M 2 , the system contains a finite number of highly populated sites, i.e., 
condensate precursors. The second stage of the non-stationary dynamics, where all but 
one of the precursors die out, also lasts a length of time of the order of the diffusive 
timescale M 2 . The whole non-stationary process of the formation of the condensate 
is therefore characterised by a single timescale r non _ st ~ M 2 . A similar scenario holds 
in the IDS geometry, the only difference being that r non _ st now scales as M 3 . Table 1 
summarises the values of the dynamical exponents z and Z, such that r non _ st ~ M z 
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and r ~ M z , where r ~ T\ oc ~ r g i is the characteristic timescale for the stationary 
motion of the condensate. The shift of the dynamical exponent by one unit in the IDS 
geometry has a common origin [7, 8]: it stems from the Gambler's ruin problem [20]. An 
analogous phenomenon is encountered for example in the coarsening law for the domain 
growth, and in the motion of a tagged particle, in ID Kawasaki dynamics [21]. 



Geometry 


z 


Z 


MF, IDAS 


2 


b 


IDS 


3 


6+1 



Table 1. Non-stationary and stationary dynamical exponents of the ZRP with 
static exponent b > 2. 



As recalled above, the non- stationary dynamical exponents are insensitive to the 
exponent b, and more generally to the statics, provided the system is in its condensed 
phase. This feature is easily understood in the context of the Markovian Ansatz 
proposed in the present work. The last stage of the formation of the condensate, i.e., 
the disappearance of the smaller of the last two precursors, implies no barrier crossing. 
In terms of the occupation of the condensate, it indeed corresponds to the transition 
from Ni to A, where the initial occupation N± of the larger precursor was already larger 
than A/2, corresponding to the top of the potential barrier. This explains why r non _ st 
is given by the diffusive timescale, both in the framework of the Markovian Ansatz and 
in the MF and IDAS geometries. 

The prefactors appearing in (6.2) demonstrate that the dynamics is faster in the 
IDAS geometry than in the MF geometry by a factor of order 2/1.26 ~ 1.59. This 
effect can be explained at least qualitatively in terms of the mean time T mean taken 
by a particle to travel between two randomly chosen sites. In the IDAS geometry, the 
motion is ballistic, implying that time is just distance, hence T mcan = M/2. In the MF 
geometry it can be shown that T mean = M. The factor 2 between these two estimates 
provides a rough explanation of the observed ratio of characteristic times. 

The present approach can be easily generalised to other classes of ZRP. One 
situation of interest is when the critical occupation probabilities fall off more rapidly 
than any power law (so that formally b = oo). Using both the expression (A. 12) of 
the probabilities in the dip region, and the Arrhenius law (5.2) (or (B.13)) for barrier 
crossing, we predict that the characteristic time scales as 

1 Pa 

/a/2 Pa/2 2 

This estimate diverges faster than any power of M, but not exponentially fast. 
The power-law prefactor, which is missed by using the Arrhenius expression r A 
instead of the full prediction r MA , is therefore irrelevant. An explicit example is 
provided by the rate (A. 13). The critical occupation probabilities decay as the 



Tma ~ ~ ~ ~ J" 
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stretched exponential (A. 14), whereas the characteristic time has a stretched exponential 
dependence in the system size M: 



The predictions obtained in this work can also be extended to other geometries, 
besides the MF and ID cases, in analogy with the case of the non-stationary dynamics of 
formation of the condensate [22]. The MF exponents (first line of Table 1) are expected 
to hold for partially or totally asymmetric dynamics in any dimension, and for symmetric 
dynamics in dimension d > 2. Finally, a logarithmic factor is expected for symmetric 
dynamics in two dimensions: r In M ~ M b In M. 

As a last remark, let us mention that it would be interesting to see whether the 
predictions made in the present work extend to the case of driven diffusive systems 
leading to condensation. If, as currently believed, these systems can be adequately 
mapped onto ZRP in the stationary state [11], one would expect this to hold. 

We end up with a critical review of references [8, 10], where the stationary-state 
motion of the condensate had already been addressed. Reference [8] contains predictions 
both on the late stage of the formation of the condensate (non-stationary dynamics), and 
on the dynamics of the condensate in the stationary state. These predictions rely upon 
the analysis of an effective symmetric two-site model. In the stationary state this model 
leads to an inverse binomial distribution. (By means of an argument equivalent to the 
Arrhenius law, the characteristic time for a macroscopic fluctuation of the condensate 
is then found to grow exponentially with the system size.) In our language, reference [8] 
predicts that the stationary occupation probabilities are exponentially small in the dip 
region, which is in strong disagreement with the power-law scaling result (5.1), itself a 
direct consequence of the exactly known finite-size stationary probabilities (2.3). The 
effective model proposed in [8] is therefore inappropriate to describe the stationary 
ergodic motion of the condensate. (One possible explanation is that the rates defining 
the effective model are derived under an implicit assumption of weak fluctuations in 
the probability currents, which might not hold when the condensate jumps from one 
site to another.) Reference [10] is devoted to the statics of finite-size systems for a 
class of one-dimensional mass transport models leading to condensation. The authors 
also give an estimate for the characteristic time r of the dynamics of the condensate 
at stationarity. This time is argued to diverge as some power of the system size in the 
regime of anomalous fluctuations (2 < b < 3) and exponentially fast in the regime of 
normal fluctuations {b > 3). This prediction is incorrect because the required probability 
is estimated by using an expression for the condensate hump outside its range of validity. 
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Appendix A. Statics of the ZRP 

This Appendix presents a reminder of the main static properties of the ZRP, and a 
derivation of a few results relevant to section 5. 



Simple facts 

The partition function Zm,n entering (2.1) reads 

m / M 



M / M \ 

Zm,n= E I[pnJ(T, N ^ N ) > 

{Ni} i=l Vz=l / 



and obeys the recursion formula 

N 

Zm,n — EPfc Zm-i,n-h- (A.l) 

fc=0 

This ensures that the stationary occupation probabilities 

/ fe , s t = PfcZ ^-^-* (A.2) 
are normalised. We have 

JV 

-^O.AT = ^AT,0) -^l,Af — PN, Z 2t N = E PkPN-k, (A. 3) 

fc=0 

and so on. Using an integral representation of the Kronecker delta function, we obtain 

^•Wd^M". (A ' 4) 

where the generating series of the weights pk reads 
P(z) = 5>*z*. 

A;>0 

Static properties of the ZRP are therefore entirely encoded in this series. In the 
thermodynamic limit (M — > oo at fixed density N/M = p), the free energy per site, 

JF= ~ »i im T7 lnZ M,N, 

M->oo M 

can be obtained by evaluating the contour integral in (A. 4) by the saddle-point method. 
The saddle-point value z depends on the density p through the equation 

z P'(z ) 



The free energy per site and the stationary occupation probabilities read 

Pk ^ 

T = p\nz - In P(z ) , f k ,st = 



(A.5) 



p(z y 
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Rates Uk — 1 + b/k: Power-law critical behaviour 

For the choice of rate Uk — 1 + b/k, the above formulas yield 

F T(A; + 6+l) io V ; fc 6 ' 

P(z) = / A- ] — du = 2 F 1 (1, 1; 6 + 1; z), (A.6) 

where 2 Pi is the hypergeometric function. The function P(z) has a branch cut at z — 1, 
with a singular part of the form 

P sg (z)«AP(l)(l-z) & -\ 

so that P(z) is only differentiable Int(6) — 1 many times at z — 1. The following values 
are of interest: 

0—1 Sill 7T0 /» 7\ 

6 46 t A -'J 

P'(l) = P"(l) = — 

K) (6-l)(6-2)' K) (6-l)(6-2)(6-3)' 

Whenever b = n > 2 is an integer, the amplitude A diverges. The singular part of the 

generating series is of the form P sg (z) « n(— 1)™(1 — z) n_1 ln(l — z). 

For 6 > 2, the system has a continuous phase transition at a finite critical density 

P^l) J_ 
^ c P(l) 6-2' 

such that the saddle point zq reaches the singular point z — 1. This critical density 
separates a fluid phase (p < p c ) and a condensed phase (p > p c ). 

Critical density (p = p c ). The occupation probabilities 

_ p fc ^ (6-i)r(6) 

/fc,st - -^jy ^ (A.SJ 

fall off as a power-law in the thermodynamic limit. The critical free energy reads 

J- c = -lnP(l) = -ln^ T . 
The second moment of the occupation probabilities, 

^ C ~h P ^ (6-2)(6-3)' 

is convergent for 6 > 3 (regime of normal fluctuations), and divergent for 2 < 6 < 3 
(regime of anomalous fluctuations). 
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Condensed phase (p > p c ). A large and finite system in the condensed phase essentially 
consists of a uniform critical background, containing on average N c = Mp c particles, 
and of a macroscopic condensate, containing on average A = N — N c = M(p — p c ) excess 
particles with respect to the critical state. 

The occupation probabilities /fe ;St accordingly split into two main contributions [10]. 
The contribution of the critical background, corresponding to small values of the 
occupation (k <C M), is approximately given by (A. 8). The contribution of the 
condensate shows up as a hump located around k — A. The hump is a Gaussian 
whose width scales as M 1 / 2 whenever p c is finite, i.e., for b > 3, whereas it has power- 
law tails and a larger width, scaling as M l ^ b ~ l \ in the regime of anomalous fluctuations 
(2 < b < 3). The weight of the condensate probability hump is approximately 1/M, 
in accord with the picture that the system typically contains a well-defined condensate 
located on a single site at any given time. 

For the investigation of the stationary dynamics of the condensate, we need to 
estimate the occupation probabilities /fe ;St all over the probability dip region lying 
between the critical background and the condensate hump, i.e., when k >> 1 and 
A — k >> 1. To reach this goal, we first estimate the partition function Zm,n deep 
in the condensed phase, i.e., for A = M(p — p c ) ^> 1. An analogous analysis, with a 
different goal, is given in [9]. In this regime, the contour integral in (A. 4) is dominated 
by the vicinity of the singular point z = 1. We therefore set z = e~ £ , anticipating that 
the relevant e will be small. Equation (A. 4) then reads 

Z M , N = P(l) M J c exp [M ((p - Pc )e + ■■■ + Ae^ 1 + •••)]. (A.9) 

In this expression, the contour C encircles the real negative axis, whereas the first dots 
stand for higher-order regular terms (e 2 , e 3 , and so on) and the second dots stand 
for higher-order singular terms. The leading contribution comes from linearising with 
respect to the leading singular term: 

* i / . ,.\ ~ AM P(1) M [ — e A . ; . 



/c 27ri 

The contour integral reads 1/(T(1 — b) A b ). We are thus left with 

M ( b \ M 

ZM,N*{b-l)T{b)-^—} . (A.10) 

(This expression also depends on M and N through A = N—Mp c .) The last exponential 
factor is exp(— MjF c ): the free energy throughout the condensed phase is equal to the 
critical one. Note the exact correspondence with the formula (4.3) for the two-site 
system for b > 1, up to the replacement of M by 2 and of A by N. Equation (A. 2) then 
yields 

^"(i-lTOp^j- (A.11) 

This estimate holds throughout the dip region (k ^> 1 and A — k ^> 1). The occupation 
probabilities (A. 11) are proportional to the product PkPA-k- The normalisation matches 
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with (A. 8) for 1 <C k <C A. The above results are fully insensitive to the (normal or 
anomalous) character of critical fluctuations. 

The natural interpretation of the result (A. 11) is that typical configurations where 
one site contains k particles in the dip region (k 1 and A — k ^> 1) are configurations 
where the remaining A — k excess particles live on a single other site. The dip region is 
therefore entirely dominated by configurations where the excess particles are shared by 
two sites. It can be checked that configurations where the excess particles are shared by 

p > 3 sites (i.e., k + k\ H V k v ^\ = A with fc> 1 and all the ki ^> 1) have a negligible 

total weight. This weight indeed comes from integrating the term (Ae b ^ 1 ) p ^ 1 in the 
right-hand-side of (A. 9). It is therefore suppressed by a power- law factor \/^p-^-i) 
with respect to the leading result (A. 11). 

The derivation of (A. 11) can be extended to the ZRP with an arbitrary critical 
mass distribution. The occupation probabilities in the dip region are still dominated by 
configurations where the excess particles are shared by two sites. They read 

A, » ^ (A.!2) 
Rates Uk = 1 + a/k a : Stretched- exponential critical behaviour 

We now turn to a brief account of the static properties of the ZRP with hopping 
rate [11, 12] 

Uk = l + ±, (A.13) 

where a is an arbitrary exponent. The situation of interest corresponds to < a < 1. 
Equation (2.2) leads to the estimate 

Pk ~ exp (^-a ^ ~ exp (- k 1 -^ . (A. 14) 

The generating series P(z) has an essential singularity at z = 1 with an exponentially 
small discontinuity. The occupation probabilities at the critical density, fk, s t = Pfe/P(l), 
decay as a stretched exponential law. 

For a large finite system, the stationary occupation probabilities in the dip region 
are given by (A. 12), i.e., in the present situation, 

f k , st ~ exp [A 1 -" - k 1 - - (A - k) 1 - 

\ 1 — a L 

Appendix B. Biased diffusion on an interval 

This Appendix gives a brief presentation of the properties of a biased random walk on an 
interval in continuous time, with arbitrary nearest-neighbour hopping rates. The main 
emphasis is on the various definitions of characteristic times and on their relationships. 
Part of the content of this Appendix is classical and can be found in the literature [23]. 
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Master equation. Sites are labeled by k — 0, . . . , N. The hopping rates are \ k to the 
right (k — > k + 1), and Hk to the left (k — > fc — 1). The boundaries are reflecting: 
A at = /io = 0. The probability fk(t) that the random walk is at site k at time t obeys 
the master equation 

= fjik+ifk+i + Afe-i/fc-i - (//fe + Afc)/ fc (1 < A; < N - 1), 
=|ii/i-Ao/ , (B.l) 

— — — An-IJN-1 — fJ-NJN- 

at 

This equation can be recast as 

—rr — Jk-i — Jk-, 
at 

by introducing the probability current through the bond (k, k + 1): 

Jk = ^kfk — Hk+lfk+1- 

Stationary state. In the stationary state, the current Jk, s t is independent of k. The 
boundary conditions impose that it vanishes. The detailed balance condition 

Afc/fc, s t = A*fc+i/fe+i,st (B-2) 
is therefore fulfilled, implying that the system is at equilibrium. Equation (B.2) yields 

/fc,st = /o,st II ' (B-3) 

1=1 V-L 

where /o,st is fixed by normalisation. The occupation probabilities can be rewritten as 

jk,st /o,st6 k , 

in terms of the effective potential§ 

k 

e=i A ^-i 

Spectrum of characteristic times. The master equation (B.l) reads formally 

where M. is the Markov matrix of the process. The spectrum of characteristic times is 
obtained by diagonalising M. Setting f k (t) = (p k e~ Et , we obtain 

Hk+l¥k+l,a + Afe_i</?fc_i ia — (/ifc + Xk)fk,a — ~E a ip k;a . (B.4) 

The ground state E = corresponds to the stationary solution ipk,o = fk,st of (B.3). 
The characteristic times r a are the inverses of the eigenvalues E a for a — 1, . . . , N. 

§ In this paper we indifferently use words describing the occupation probabilities themselves (hump, 
dip), and the associated effective potential (valley, barrier). 
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Mean first-passage time. Let T kj0 be the mean first-passage time by the origin of 
the particle with initial position k, i.e., the mean time it takes the particle to reach 
site 0. This mean first-passage time obeys the following backward equation, obtained 
by conditioning on the first step of the walker during the infinitesimal duration dt, 

T kfl = dt + (1 - (A fc + ii k )dt) T kfi + \ k dtT k+lfl + fi k dtT k _ lfi , 

hence 

(A*, + /ifc)T fci0 — \ k T k+ ifl — /ifcT fc _ 10 = 1. (B.5) 

The backward equation (B.5), with boundary condition T 0j o = 0, can be solved 
in closed form. Indeed the differences d k = T kt0 — T k _i t0 , defined for k > 1, obey 
the two-term equations fi k d k — \ k d k+ i = 1. The homogeneous equations, without 
the right-hand-side, have the solution d k = l/(fJ>kfk,st)- Looking for a solution of the 
inhomogeneous equations in the form d k = c k / (fJ,kfk,st), we obtain the difference equation 
c fc - c fc+ i = f k , s t, with formally c N+1 = 0, hence c k = J2eL k fe,st, and finally 

k y N 

T k fl = —7 fmjBt- ( B -6) 

e=i MJ^st m=e 

We now define the following two crossing times: T L = T Nj0 (resp. T R = T 0jN ) is the 
mean time it takes the particle to cross the system from right to left (resp. from left 
to right), i.e., the mean time taken by the particle at initial position N to reach site 
(resp. at initial position to reach site N). These crossing times read 

N -y N N ^ t-l 
Th — ^ 7 X! f™,st, — 7 X! fm,st- (B-7) 

£ =1 Heje,st m= £ £ = i fi>ej e,st m= o 

The sum of these two times is the return crossing time R, that is the mean time taken 
by the particle to cross the system from one end of the system to the other and back: 

N i JV-1 I 

T L + T R = £> = £— -= . (B.8) 

e=1 n>eje,st £ =0 Mjt,st 

The expressions derived so far are quite general. They can be simplified in the 
following situations of interest. 

The symmetric random walk. The case of the usual symmetric random walk corres- 
ponds to uniform unbiased hopping rates = fi k = A at all sites k. Equation (B.4) 
becomes the discrete Laplace equation 

A(y?fc+i + (pk-i - 2ip k ) = -Eip k , 

with boundary conditions </?_i = ip and ip^ — <Pn+i- The stationary state is uniform: 

f 1 

Jk,st 



N+ 1 

The other eigenvalues of the Markov matrix read E a = 4A sin 2 [a7r/(2(iV + 1))]. 
Using (B.6), the mean first-passage by the origin reads T k $ = k(2N + 1 — fc)/(2A), 
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hence the crossing times have the simple expression Tl = Tr = (N(N + 1))/(2A). In 
the limit of a large interval (N ^> 1), we have therefore 

N 2 N 2 
Tl = Tr -2A' 

hence 

T L = T R ^yr 1 . (B.9) 

The crossing times and the main characteristic time grow proportionally to the diffusive 
timescale N 2 . 

The case of a high potential barrier. The situation of most relevance in the context of 
this paper is when the effective potential V k exhibits a high barrier, i.e., a pronounced 
maximum Vb I at some position k = B. The barrier defines two potential valleys, the 
left one (0 < k < B), and the right one (B < k < N). For simplicity, we assume that 
the effective potential is small everywhere except in the region of the barrier (k « B). 
The respective stationary weights of the valleys read 

B-l N 

Ph — /fc,st, -Pr — X/ /fc,st- 

fc=0 k=B+l 

These two weights sum up approximately to unity since fs,st ~ g~ Vb is negligible. 
The expressions (B.7) for the crossing times simplify in this regime. Sums over i are 
dominated by the contributions of sites t ~ B close to the top of the potential barrier. 
We thus obtain the estimates 

Tl ~ PrR, T r ^P l R, (B.IO) 

since the return crossing time R (see (B.8)) is also dominated by the vicinity of the top 
of the potential barrier. 

The main characteristic time T\ = 1/Ei is also expected to be very large in the 
presence of a high potential barrier. It can be estimated as follows. Let ipk t i be the 
corresponding eigenvector. Introducing the ratios 

n-k — -j. — , 
Equations (B.4) are equivalent to 

Jk — Jk-i — Exfk^hk, h k — h k+ i = . (B.I I) 

fJ'k+l Jfc+l,st 

First, the current J k can be neglected, so that h k is constant. This scheme holds as long 
as the /fc jSt are not too small, i.e., separately in each valley. We have therefore 

h k ^a L (0 < k < B), h k ^a R (B < k < N). 

The sum rule Y.k=o i Pk,i = then implies ~ Pr, Or ~ —Pl, up to an overall 
multiplicative constant. Second, the first equation in (B.ll) then yields the current, 
to first order in E\. 

J k » P K E, 2 fe >Bt (0 < k < B), J k ^ J2 kst (B <k < N). 

e=o t=k+i 
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Both expressions consistently give Jb ~ P r Pl-&l- Third, using the second equation 
of (B.ll), we have 

1 ps P L + P R ps a L — a R ps /i — h N ps J B P ps PrP^ExR. 
We thus finally obtain the following estimate for the main characteristic time: 

t^PrPlR, (B.12) 
as well as the relationships 

n » P L T L ps p r t r , I » i- + 2- 

Tl -?L Tr 

between the main characteristic time and the crossing times. 

If the potential barrier is very high, i.e., if f min is very small, a crude estimate of 
the characteristic times can be obtained by means of a saddle-point evaluation of the 
sums in (B.7): 

n ~T L ~T R ~ ~e Vs . (B.13) 

J min 

The Arrhenius law (5.2) is thus recovered. The corresponding activation energy is just 
the top of the potential barrier, Vb — — In / min . 

The full expressions (B.10) and (B.12) however provide more accurate predictions 
for the characteristic times in every situation pertaining to the high-barrier regime. 
They indeed lead to the correct absolute prefactor to the estimate (B.13), which usually 
involves a power of N. For instance it is used in the present work in situations where the 
characteristic times grows as a power of N. The above mentioned prefactor is needed 
in general in order to correctly predict this power. Finally, let us mention the following 
self-consistent viewpoint on the relevance of the high-barrier regime: it is expected to 
hold whenever the predicted characteristic times grow much faster than the diffusive 
timescale N 2 . This viewpoint is illustrated by the case of two sites (see section 4): the 
high-barrier regime corresponds to b > 1. 

For the case of a symmetrically biased diffusion problem, the two valleys are equally 
weighted (Pl — Pr — 1/2). The crossing time T = T R = Tl and the main characteristic 
time Ti are thus given by 

n « | = f ■ (B.14) 

For a very asymmetric problem, with a small right valley of weight Pr <^ 1 and a large 
left valley of weight P L ~ 1, the crossing time to the large valley is close to the main 
characteristic time (Tl ~ n), whereas the crossing time to the small valley is much 
longer, as it reads approximately T R ps t^/Pr, in general. 
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